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ABSTRACT 

The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic 
Cap, using the QOPrime wide-held imager on the 2.3 Bok telescope at Kitt Peak. The survey observations started in 
2010 and ended in 2013. The hnal survey area is about 5000 deg^ with a median 5cr point source limiting magnitude of 
^ 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric 
and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on 
objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread 
function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model 
magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS 
photometric catalog with those of the SDSS and CanadaFranceHawaii Telescope Legacy surveys. 

Subject headings: surveys — catalogs — techniques: image processing — techniques: photometric 


1. introduction 


^he South Galactic Cap u-band Sky Survey (SCUSS; see 
J^r^hou et al. 2015, in preparation for more detail) is a 
imaging survey in the South Galactic Cap in it-band 
an effective wavelength of 3538 A. It is an international 
Serative project between the National Astronomical Ob- 
sm^tories of China and the Steward Observatory of the 
I'ersity of Arizona. The survey utilizes the 2.3m Bok 
scope at Kitt Peak. The camera, installed at the prime 
f^ojjS, provides a field of view (FOV) of about 1 deg^. The 
adjured filter is similar to the u band of the Sloan Digital 
Sl^ Survey (SDSS; York et al. 2000). The SCUSS project 
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started in the summer of 2009 and began its observations 
in fall 2010. Survey observations were completed in the fall 
of 2013. The final survey area is about 5000 deg^ with uni¬ 
form imaging depth, of which more than 75% is covered by 
the SDSS footprint. The median imaging depth for point 
sources is about 23.2 mag at a signal-to-noise ratio (S/N) of 
5 with a 5 minute exposure time. Table 1 gives the basic 
characteristics of SCUSS. 

Table 1: Survey Summary 


Telescope 

2.3 m Bok telescope 

Site 

Kitt Peak in Arizona 

CCD 

2x2 4kx4k CCD array 

Field of view 

l°.08xl°.03 

Filter 

u (3538A) 

Integral time 

300 s 

Magnitude limit 

23.2 mag 

Survey area 

~5000 deg2 

Observation period 

2010-2013 


The main goal of the survey is to supply input photomet¬ 
ric catalogs to select spectroscopic targets for the Large Sky 
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Area Multi-Object Fiber Spectroscopy Telescope (Cui et al. 
2012). In addition, by combining with other bands in large- 
scale photometric surveys, such as the SDSS and Panoramic 
Survey Telescope & Rapid Response System (Pan-STARRS 
Kaiser 2004), the survey data can be used for a wide range of 
scientific investigations, such as Galactic structure, Galactic 
extinction, galaxy photometric redshift, galaxy star forma¬ 
tion rates, and stellar populations of nearby galaxies. 

This paper describes the data reduction pipeline specially 
designed for the SCUSS survey. There are some instrumen¬ 
tal issues that need to be specially handled, such as issues 
in overscan, substructures in bias, and crosstalk. In addi¬ 
tion to detecting sources on SCUSS images, we also provide 
photometry for SDSS objects using consistent object param¬ 
eters. Section 2 introduces the SCUSS survey and related 
facilities. Section 3 describes the basic image processing. 
Astrometric and photometric calibrations are presented in 
Sections 4 and 5, respectively. Section 6 provides image- 
quality statistics of the observations. Section 7 describes 
image stacking. Sections 8 and 9 present photometry meth¬ 
ods and a comparison with other surveys, respectively. A 
summary is given in Section 10. 

2. The survey and facilities 

SCUSS is a u-band imaging survey in the northern part 
of the southern Galactic Cap. The survey originally covered 
the region of Galactic latitude b < —30° and equatorial lat¬ 
itude 6 > —10°, with a total area of about 3700 deg^. It 
was further extended to the Galactic anti-center region and 
the extra area covered by the SDSS, with a final survey area 
of about 5000 deg^ (X. Zhou et al. 2015, in preparation). 
Normally, there are two exposures for each field and the to¬ 
tal exposure time is 5 minutes, which generates images 1-1.5 
mag deeper than the SDSS u band. In the following, we give 
a general description of the telescope, camera, and filter that 
were used in SCUSS. 

2.1. Telescope 

The Bok telescope ^ is a 90 inch (2.3 m) telescope op¬ 
erated by the Steward Observatory of the University of 
Arizona. It is located on Kitt Peak, whose latitude is 
-|-30°57'46".5 and longitude is 111°36'01".6W. The eleva¬ 
tion is about 2071 m and the typical seeing is about 1".5. 


^http://James.as.arizona.edu/ psmith/90inch/90inch.html 


The telescope runs year-round except on Christmas Eve and 
during a maintenance period in August. The accuracy of 
the absolute pointing of the telescope is recorded as 3" over 
the entire sky. 

2.2. Camera 

An imaging system, named 90Prime, is deployed at the 
prime focus (corrected focal ratio: f/2.98; corrected focal 
length: 6829.2 mm). The detector is a CCD array con¬ 
sisting of four 4kx4k backside-illuminated CCDs. They are 
STA2900 CCDs made by Semiconductor Technology Asso¬ 
ciates, Inc. and backside processed at the University of Ari¬ 
zona Imaging Technology Laboratory. These CCDs have 
been optimized for the u-band response, giving a quantum 
efficiency at u band of about 80%. Figure 1 displays the 
layout of the CCDs on the focal plane. The edge-to-edge 
FOV is about 1°.08x1°.03. The pixel scale is 0".454. There 
are inter-CCD spacings in the center of the array: 166" in 
right ascension and 54" in declination. Each CCD is read 
out by 4 amplifiers located at the corners. Each detector 
has 4096x4032 physical pixels and 20-row pixels of overscan 
for each amplifier. The full well is about 90,000 electrons or 
65,000 data numbers (DNs). The current system gain is set 
to be about 1.5 electrons per DN. The dark current is about 
7 electrons per pixel per hr. The readout time is about 30 
s and the average readout noise is about 8.8 electrons. The 
response non-uniformity with the u filter, which is the stan¬ 
dard deviation divided by the mean of a u band flat field 
image, is about 2%. Table 2 summarizes these parameters 
of the detector. In 2010, CCD #4 had problems and only a 
quarter of this CCD could be used. CCD #2 and CCD #3 
had relatively large readout noise. The camera was upgraded 
in 2011, at which time these two CCDs were replaced with 
new detectors. CCD ^ 4 was replaced and swapped with 
CCD #1. New video preamps were added to all four CCDs 
in order to reduce crosstalk and improve noise immunity. 

2.3. Filter 

The SCUSS u filter is similar to the SDSS u band. Figure 
2 displays both SCUSS and SDSS system response functions. 
The SCUSS u response curve includes the Hlter transmission, 
the CCD quantum efficiency, and the atmospheric extinction 
at the typical airmass of 1.3. The adopted atmospheric ex¬ 
tinction is the same as that of the SDSS, which is based 
on the standard Palomar monochromatic extinction coefh- 
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Fig. 1.— Layout of the CCD array. There are four CCDs: 
CCD #1, #2, #3, and #4. Each CCD has four amplifiers 
highlighted by dashed lines. Overscan regions are also indi¬ 
cated. 


Table 2: Detector Parameters 


CCD number 

4 (in 2x2 array) 

CCD gaps 

166" in R.A. and 54" in deck 

CCD size 

4096x4032 (extra 40 rows of overscan) 

Full well 

90,000 electrons/pixel (65,000 ADUs/pixel) 

Amplifier number 

4 

Quantum efficiency 

~80% at 35 OOA 

Average dark current 

7.2 electrons/pixel/hr 

Response non-uniformity 

2% at u band 

Average readout noise 

8.8 electrons 

Readout time 

~30 s 

Average gain 

1.5 electrons/DN 

Plate Scale 

0".454 (15 Atm) 


effective wavelength and bandwidth of the SCUSS u band are 
about 3538 A and 345 A, respectively. The FWHM is about 
520 A. The effective wavelength, bandwidth, and FWFIM of 
the SDSS u filter is 3562, 385, and 575 A, respectively. The 
SCUSS u band is slightly bluer than the SDSS filter. In the 
rest of this paper, we will use the symbol of u* to refer to the 
SCUSS u band and the symbols of u, g, r, i, and z for the 
five SDSS bands. Throughout the paper, objects are clas¬ 
sified as point-like or extended throughout this paper using 
the SDSS star-galaxy separation unless otherwise specified. 



Fig. 2.— Response curves of the SCUSS u* and SDSS u 
filters. Both curves include the atmospheric extinction at 
an airmass of 1.3 and normalized to their maximums. 


■ 3 . Basic image processing 


cients but here scaled to the elevation of Kitt Peak assum¬ 
ing an exponential scale height of the atmosphere of 7000 m 
(Doi et al. 2010). The effective wavelength and bandwidth 
of the filter are respectively defined as AeS = 

Aeff = ’ '"^here i?(A) is the filter response curve. The 


3.1. Image division and overscan correction 

OOPrime has four CCDs and each CCD is read out by 
four amplifiers. Thus, every exposure file includes 16 FITS 
extensions. We split the file into four smaller FITS images, 
each of which represents one of the four CCDs. Overscan 
lines (40 pixels) are moved to the right side of the image. 
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We compute the median of the overscan columns for each 
amplifier and subtract it from the raw frame. There are 
some subtle issues about the overscan that require special 
care. For example, when a bright star is located right be¬ 
side the overscan region, nearby overscan pixels are contam¬ 
inated. Occasionally, some brighter stripes appear within or 
close to either end of the overscan region. The overscan in 
these regions needs to be interpolated or extrapolated. After 
overscan subtraction, frames are rotated counter-clockwise 
by 90° to keep north up and east left. These frames are 
then trimmed to the size of 4096x4032 pixels. 

3.2. Dark Current 

The average dark current of the 90Prime imager is about 
7.2e per hr. For an exposure of 300 s, it is ignorable (about 
0.6e) relative to the readout noise of about 8.8e, so we did 
not take dark exposures. 

3.3. Bias 

Bias frames are taken before and after the scientific ex¬ 
posures. A total of 20 bias frames are obtained each night. 
After combining them, there are large-scale structures, espe¬ 
cially for CCD ^2 and #3 as shown in Figure 3. Counts in 
some peaks of the structures can be more than 10 DN. These 
structures affect the accuracy of photometry significantly, es¬ 
pecially in u band, because counts of scientific frames in this 
band are low. The bias structures are stable for a long time 
(month-to-month variation is about 0.25 ADU), so we de¬ 
rive a median ” super” bias by combining all the bias frames 
taken within a month. All overscan-subtracted frames are 
corrected by this ’’super” bias. 

3.4. Flat-fielding 

Flat-fielding removes the instrumental signatures from 
raw frames using several exposures taken with the telescope 
facing a uniform light source. At the beginning and end of 
each night, a total of 20 dome flats are taken with the dome 
closed and the telescope pointing to a white screen. This 
screen is illuminated by UV lamps of Philips MasterColor 
Ceramic Metal Halide ED-17. Usually, a 6 s exposure gen¬ 
erates about 20,000 DN on the CCD. Although dome flats 
have extremely high S/Ns, there are some disadvantages to 
applying them in flat-fielding. First, lamps are point sources 
and the scattered light on the screen might be uneven in the 
radial direction. Second, the optical path of the light from 



Fig. 3.— ’’Super” bias image of four CCDs in 2011 Decem¬ 
ber. The layout of these CCDs is the same as that shown in 
Figure 1. 

the screen is not the same as the light from night sky. In 
addition, we also find that the gain of each amplifier changes 
considerably during the whole night. When using only dome 
flats, we find that average sky backgrounds in the four am¬ 
plifiers of each CCD can be quite different. 

Twilight flats are an alternative way and regarded as be¬ 
ing more uniform than dome flats. They are obtained by 
observing the sky during evening and morning twilight. The 
brightness of the twilight sky changes rapidly and it strongly 
depends on the weather. It is difficult to get enough good 
high S/N flats. The problem of gain variation also persists. 

Raw science images also contain the flat-held character¬ 
istics. By combining all the science frames with outliers 
rejected, we can obtain a ’’super” sky hat. This kind of hat 
has more advantages than the two types of hats above. The 
light of night sky ’’super” hat comes through the same op¬ 
tical path as the observed objects. The CCD gain is the 
average of all science frames and it is closer to the real-time 
variation during the night. There are some structures in the 


4 



’’super” flat that present a different sensitivity amplitude 
relative to the dome flats possibly due to a different light 
path. About 150 science images on average are observed 
each night. The average level of the sky background in each 
image is about 150 ADU. As a result, the ’’super” sky flat has 
an average count of about 23,000 ADU. We present the sky 
flats on 2011 December 28 as an example in Figure 4. After 
applying either the dome flat (higher S/N) or super sky flat, 
we find that there is less than a 0.2% of difference in the sky 
background fluctuations. The resulting error due to lower 
S/N of the ’’super” flat is negligible relative to the large sky 
background fluctuation of about 12.7 ADU and CCD read¬ 
out noise of 8.8 electrons. Thus, we use the ’’super” sky flat 
to do flat-fielding corrections of the science frames. 



Fig. 4.— ’’Super” sky flat image of four CCDs on 2011 
December 28. The layout of these CCDs is the same as that 
shown in Figure 1. 


3.5. Crosstalk 

Crosstalk frequently occurs in multi-channel CCD chip 
read-out. When there are saturated objects in one of the 
CCD amplifiers, it can be best seen as mirror images in 
other amplifiers. The crosstalk usually causes contamina¬ 
tions across the output amplifiers at the level of 1:10,000 
(Freyhammer et al. 2001). The output counts of each ampli¬ 


fier are the sum of the true counts from the sky and a small 
fraction of counts from other amplifiers. The crosstalk signal 
can be either positive or negative and should be corrected 
for high-precision astronomical photometry. 

The OOPrime camera has four CCDs and each CCD has 
four amplifiers. The crosstalk effect is clearly seen in our 
SCUSS raw images with saturated stars. The crosstalk sig¬ 
nal is positive. We characterize the effect assuming that it is 
additive and proportional to the number counts of other am¬ 
plifiers. The proportionality coefficients are relatively stable 
in SCUSS images. We correct crosstalk using the following 
prescription: a series of images with bright stars appearing in 
one of the CCD amplifiers is selected; the proportionality co¬ 
efficients are estimated by comparing the mirror signals with 
their original signals; the crosstalk signals in one quadrant 
are removed with the corresponding coefficients of the other 
three quadrants. The average coefficient for the SCUSS im¬ 
ages is about 2:10,000. Figure 5 illustrates the crosstalk and 
the performance of its correction. The crosstalk signals are 
clearly seen as mirror images (arrows) of a bright star that 
is marked with a circle. 

4. Astrometric calibration 
4.1. Astrometry by UCAC4 

Astrometric solutions are derived by cross-identifying ob¬ 
jects in science frames with the Fourth US Naval Obser¬ 
vatory CCD Astrograph Catalog (Zacharias et al. 2013, 
UCAC4). UCAC4 is an all-sky star catalog that is com¬ 
plete in R band down to 16 mag. It contains proper motions 
for most stars. An approximate first-order guess of the as¬ 
trometric solution for each CCD is first estimated by using 
UCAC4 objects in the same area with positions revised ac¬ 
cording to their proper motions. Then, the projection center 
of the telescope can be calculated using this rough solution. 
A more accurate astrometric solution is then derived with 
the projection center and a radial second-order correction 
for the focal plane distortion. This distortion term is the 
same for all frames in our survey. 

The left panel in Figure 6 shows the distribution of the 
astrometric errors in the plane of R.A. and deck differences 
between SCUSS and UCAC4. It includes objects in one 
exposure of a randomly selected field. About 1800 objects 
are matched with UCAC4. The average position offsets are 
0".002 and 0".004 for R.A and deck, respectively. The Icr po- 
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Fig. 5.— Left: a SCUSS raw CCD image with visible crosstalk signals. The arrows point to the mirror images of the bright star 
marked with a circle. Right: the image after correcting the crosstalk. 


sition errors in R.A. and deck are about 0".128 and 0".120. 
The right plot of Figure 6 illustrates the distortion of the 
focal plane. The plate scale varies from the center (0".455) 
to the edge (0".447) of the FOV. The pixel area near the 
corner is about 3.3% smaller than that of the central pixel. 

4.2. External astrometric errors 

External astrometric errors are estimated by matching 
SCUSS objects with the UCAC4 catalog. About 165 ob¬ 
jects on average for each CCD are crossing-identified and 
the global external position error is about 0".13±0.02. The 
mean external astrometric offsets and RMS errors for R.A. 
and deck are 

= -0".0012 ± 0.0079, M = -0."0015 ± 0.0071, 

CTAa = O'Mlll ± 0.0178, fJAi = 0".1094 ± 0.0182, 

where cr corresponds to the 68.3% confidence level. 

Figure 7 displays the external coordinate offset and RMS 
error as functions of RA and deck The offsets and RMS 
errors are uniform over the entire survey. 

4.3. Internal astrometric errors 

Internal astrometric calibration errors are estimated by 
using the cross-identifications of UCAC4 sources inside the 
overlapping areas between two adjacent exposures. Each 
field is observed with two exposures dithered by half a CCD 
size, so every CCD frame is covered by four other frames, 
whose astrometric solutions are independently derived. We 
calculate the internal astrometric uncertainties using the ob¬ 
jects in common to both frames. The global average internal 


error is about 0".09±0.03. The mean internal R.A. and deck 
offsets and RMS errors over the whole survey are 

Aa = 0".0009 ± 0.0102, AS = 0."0004 ± 0.0113, 

CTAa = 0".0728 ± 0.0295, cta^ = 0".0707 ± 0.0269. 

5. Photometric Calibration 
5.1. External calibration 

The SDSS imaging survey has covered more than one 
third of the sky within both northern and southern Galactic 
caps. It provides photometric catalogs of about 5200 square 
degrees in the SGC. More than 21 million objects in this area 
are recorded. The exposure time for each band is about 54 
s. The u magnitude limit with 95% completeness for point 
sources is about 22.0 mag. The photometric calibration ac¬ 
curacy in u band is about 1.3%. The ninth data release (Ahn 
et ak 2012, DR9) is utilized to make photometric calibrations 
of the SCUSS images. 

Because of gain and observational condition variations, 
individual zeropoints must be determined separately for 
each frame. Aperture photometry is performed by using 
DAOPHOT (Stetson 1987) on the bias-subtracted and flat- 
fielded images. We choose an aperture radius of about 
7".3 (16 pixels), similar to that of 7".43 adopted by the 
SDSS photometric calibration. The aperture diameter is 
about 7 times the typical seeing FWHM, which is large 
enough not to be affected by aperture correction. The trans¬ 
formation equation between the instrumental magnitudes 
^iinst = —2.51ogioADU to the SDSS calibrated magnitude 
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Fig. 6.— Left: distribution of astrometric errors in the plane of R.A. and decl. differences between SCUSS and UCAC4 positions 
for a randomly selected field. The red cross denotes the coordinate origin. Histograms showing the marginal distributions of the 
R.A. and decl. errors are plotted in the upper and right panels, where the red curves show the Gaussian fits. The means and 
standard deviations of the two distributions are also presented in the plot. Right: pixel scale variation in the detector plane. 
The layout of the CCD array (dashed boxes) is also plotted. Numbers labeling the contours give the pixel scales in arcseconds, 
which are also indicated by the colored surface. 
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Fig. 7.— SCUSS R.A. and decl. offsets and RMS errors relative to UCAC4 as functions of the equatorial longitude and latitude. 
The two curves in each panel show ilcr scatter around the average along the x axis. The R.A. is transformed to the range of 
(- 180 °, 180 °)- 


is simply given by 

u = Uinst + c + kX + f{u - g) + f'{u - g, X), (1) 

where c is the instrumental zeropoint, k is the atmospheric 
extinction coefficient, X is the airmass, u — 5 is the SDSS 
color, and f{u — g) is the color term. The final term of 
f'{u — g,X) is related to the color effect due to the atmo¬ 
spheric extinction. The central wavelength of SCUSS u* 
band at an airmass of 1.0, 1.3, and 1.5 is estimated to be 
3534 A, 3538A, and 3542A. The photometric effect in the 
color terms for most main-sequence stars at different air- 
masses is less than 0.3%, estimated using theoretical stel¬ 
lar spectral libraries. So, we ignore the second-order color 
term. Since the standard stars are the common SDSS ob¬ 
jects within the same area of each image, the term solely 
related to the airmass is constant. Thus, Equation (1) can 
be simplified as 

U* = U- f{u - g) = Minst + C, (2) 

where u is the calibrated SCUSS u-band magnitude and C 
is still termed as the ’’photometric zeropoint”. The color 
term is estimated as follows: (1) SCUSS instrumental mag¬ 
nitudes are calibrated using the photometric magnitudes of 


common SDSS objects without considering the color term; 

(2) magnitude differences between these calibrated SCUSS 
magnitudes and SDSS magnitudes are calculated as a func¬ 
tion of the SDSS u — g color; (3) an approximate system 
transformation formula is fitted by a two-order polynomial: 

u*-u = -f{u-g) 

= 0.1459- 0.1957(m- 5 ) 

-b0.0591(u-g)4 (3) 

for 0.8 < u — g < 2.7 (see Figure 8 ). This color term is only 
used for the photometric calibration. Here the SDSS is used 
for the zeropoint only, and the SCUSS photometry on its 
native system is defined to match the SDSS a.t u — g = 0. We 
choose stars with 16.0 < u < 20.5 to eliminate objects that 
are saturated or have large photometric errors. Equation 

(3) is then applied to transform u to SCUSS u*. Following 
Equation (2), we measure the final zeropoint C by iteratively 
rejecting outliers. 

The CCD gain varies slightly during the observation. Fur¬ 
thermore, CCD gains of the four amplifiers did not change 
synchronously. In addition, there is a photometric response 
non-uniformity on the flat-helded images. It is possibly 
caused by the focal plane distortion and scattered light re- 











Fig. 8.— Magnitude difference between the SCUSS u* and 
SDSS It as a function of the SDSS u — g color. The points 
are for point-like objects with 16 < u < 20.5. The error bars 
show the averages and Icr errors within different color bins. 
The red curve is the fitted transformation equation as shown 
in Equation (3). 
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Fig. 9.— Non-uniformity pattern used to correct the flat- 
fielded and calibrated images. The layout of these CCDs are 
the same as that shown in Figure 1. 



fleeted in the optical system (Regnault et al. 2009; Betoule 
et al. 2013).Thus, we perform photometric calibration for 
each amplifier of each CCD. There are about 47 stars in 
each amplifier on average to derive the zeropoint, whose ac¬ 
curacy is estimated to be about 0.01 mag. The zeropoint 
here is expressed in units of mag for lADU s“^. After de¬ 
riving independent zeropoints for each amplifier, we find a 
small residual pattern in the sensitivity as derived from ob¬ 
jects observed in different positions in the field as shown in 
Figure 9. We use this residual map to further refine the 
flat-held. 

5.2. Photometric response differences of the four 
amplifiers 

The average zeropoints of four amplifiers are shown in 
the left panel of Figure 10 as a function of time. Usually 
during the night, we scan the sky from east to zenith and 
then to west. The airmass and corresponding atmospheric 
extinction first decreases and then increases. Therefore, ze¬ 
ropoints present nightly variations. From this figure, it can 
be seen that the weather conditions over 2010 are worst. 
The right panel of Figure 10 shows the relative zeropoint 
variation with time for each amplifier of CCD ^1, compared 


with the average zeropoint of all amplifiers. The relative 
zeropoint variations do not keep the same value and some¬ 
times two amplifiers present opposite variations. The detec¬ 
tors have bad qualities in 2010 so that the relative zeropoint 
changes more notably with time. Thus, it is critical that we 
obtain the photometric solutions independently for the four 
amplifiers. 

Table 3 gives the zeropoint scatter of four amplifiers 
for each year and each CCD. The zeropoint differences 
among amplifiers originate from gain variation, photomet¬ 
ric response non-uniformity as mentioned before, and the 
photometric calibration error. The overall zeropoint scat¬ 
ter of four amplifiers around the averages are 0.012±0.006, 
0.016±0.007, 0.015±0.009, and 0.011±0.006 for CCD #1, 
#2, ^3, and #4, respectively. The general response differ¬ 
ences of four amplifiers in all CCDs are less than 1.5% except 
the CCD ^3 of 2010, which is about 2.8%. 

5.3. Internal calibration 

Each SCUSS held is observed twice and overlaps with the 
adjacent fields. Thus, each CCD image can be calibrated 
using the common objects cross-identified in surrounding 
images. The calibrating process begins with balancing the 
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Fig. 10.— Left: average zeropoints of the four amplifiers for each CCD as a function of time. The data points correspond to the 
CCD images, which are sorted by the observation time. The zeropoints are in units of mag for 1 ADU s“^. The data for CCD 
#1, #2, and #3 are shifted by 2.0. The dashed lines separate the years of observations. Right: relative zeropoints compared 
with the above average for each amplifier of CCD ^1 as a function of time. The data for three amplifiers , Amp. ^1, #2, and 
#3, are shifted by 0.2. 
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Table 3: Average zeropoint scatter of four amplifiers for each 


CCD 


CCD 

year 

#1 

#2 

#3 

#4 

2010 

0.011±0.007 

0.015±0.009 

0.028±0.009 


2011 

0.013±0.006 

0.015±0.006 

O.OllitO.006 

O.OlOitO.006 

2012 

0.013±0.006 

0.013±0.006 

O.OlOitO.006 

O.OllitO.006 

2013 

O.OlOitO.006 

0.020±0.007 

0.017±0.006 

O.OllitO.006 

All 

0.012±0.006 

0.016±0.007 

0.015±0.009 

O.OllitO.006 


zeropoints of images covered by the SDSS and then transfer¬ 
ring the photometric solutions to the images out of the SDSS 
footprints. The calibration iterates until the whole grid of 
photometric solutions finally converges. Figure 11 shows the 
magnitude difference distribution between measurements for 
objects that are observed twice. The histograms in black 
and red are the distributions of the same objects with the 
SCUSS it-band magnitude calibrated externally and inter¬ 
nally, respectively. The dispersion by external calibrations 
is about 0.028 mag, while that by internal calibrations is 
about 0.025. It seems the internal calibration is at least as 
good as or even a little better than the external calibration. 

6. Image Quality Statistics 

Both weather and instrumental status affect the image 
quality. The camera was regularly updated after each ob¬ 
servation season, so its performance was improved gradu¬ 
ally. The weather and night sky conditions strongly affect 
the imaging depth of the survey. We can measure the char¬ 
acteristics tracing the image quality, such as airmass, sky 
background brightness, seeing, and photometric zeropoints. 

6.1. Seeing 

The seeing is estimated by the FWHM measurements of 
isolated and bright point sources. The seeing distribution 
is presented in Figure 12a. The best seeing is 1".2 and the 
overall median seeing is about 2".0. The Bok telescope is 
located in the trough between two peaks of the mountain, 
so the wind speed is usually larger than other places on the 
same site, which has an effect on the seeing. 

6.2. Airmass 

Figure 12b shows the airmass distribution of all survey 
images. The median airmass is about 1.28, and this is used 
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Fig. 11.— Magnitude difference distributions between two 
measurements for objects observed more than twice. The 
black histogram is externally calibrated by the SDSS cata¬ 
log, while the red one is internally calibrated. Only bright 
objects with 16 < u <19 are selected. 

to determine the typical u-band filter response as shown in 
Section 2.3. 

6.3. Photometric zeropoint 

The variation of the photometric zeropoint mainly reflects 
the change in atmospheric transparency. The distribution 
of the photometric zeropoints is plotted in Figure 12c. The 
mean photometric zeropoint is 23.81 mag for lADU s“^. 
The median zeropoints for airmasses of 1.0, 1.2, and 1.4 are 
23.93, 23.85, and 23.75, respectively. There are some images 
with low zeropoints, most of which were observed when the 
weather was cloudy. 

6.4. Sky brightness 

The sky background brightness is an important parameter 
to quantify a ground-based observation station. Compared 
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with other light sources, the artificial light pollution from 
nearby cities is more serious at Kitt Peak. The distribution 
of the it-band sky brightness at Kitt peak is shown in Fig¬ 
ure 12d. There are some images taken when the moon is 
above the horizon. The average night sky brightness of all 
observations is about 22.05 mag arcsec“^. The moonless me¬ 
dian sky background at zenith is about 22.37 mag arcsec“^, 
which is comparable to that of the Apache Point site (22.1 
mag arcsec”^). Note that the calibrated sky brightness uses 
the frame zeropoints that implicitly include atmospheric ex¬ 
tinction, so the actual sky brightness is darker because the 
u-band atmospheric extinction coefficient is about 0.5, as 
estimated with the observations taken during photometric 
nights. 

7. Image Stacking 

7.1. Quality Control 

Most SCUSS fields are observed under good weather and 
moderate seeing conditions. There are some cases with bad 
image quality: (1) high sky background due to observations 
during astronomical twilight or when moon was up; (2) large 
seeing due to strong wind; (3) low atmospheric transparency 
due to cirrus in the FOV; (4) bad focus of the CCD camera. 
In 2013, we spent one and half observation runs re-observing 
most bad-quality fields. 

To ensure the homogeneity of the imaging depth, and 
completeness of the SCUSS survey, we keep only the images 
with seeing <3".0, sky background ADU <500, and photo¬ 
metric zeropoint >22.56 mag. About 92.6% of the survey 
area is covered by these images. For the remaining area, 
we take images with seeing <3".0 (3.6% area), regardless of 
the sky brightness and photometric zeropoint. If none are 
available, we take all remaining images, which covers about 
3.8% of the total area. 

7.2. Image Resampling and Stacking 

Based on the central coordinates of each field, we stack 
related single-epoch images to form a combined image. We 
first project and resample the single-epoch images to a grid 
with a fixed pixel scale of 0".454. The grid has 8640x8200 
pixels, covering a sky area of l°.090x 1°.034. If a single¬ 
epoch images contributes less than less than 128x128, it is 
not included in the stacking process. For each pixel of the 
grid, there are more than four related pixels in each single¬ 


epoch images. These pixels are regarded as having the same 
size after being fiat-fielded. We calculate the fractions of 
these pixel areas that are covered by the grid pixel and sum 
them to conserve flux. 

We subtract the sky backgrounds from the resampled im¬ 
ages. Their photometric zeropoints are converted to linear 
flux weights. The remaining signal after removing back¬ 
grounds are weighted by these weights and then co-added. 
If there are more than three pixels involved in stacking, a 
cosmic-ray rejection is implemented using a sigma clipping 
algorithm. We redo the flux calibration for each stacked 
image with the SDSS DR9 catalog to derive the final pho¬ 
tometric zeropoint. This zeropoint is approximately 29 mag 
for 1 DN. In addition, a mask image for each field is also 
generated. Each mask pixel presents the number of images 
that are actually involved in the flux co-adding. The mask 
value is reduced by one for each epoch in which a pixel is 
bad, saturated, or blank. Figure 13 shows a typical stacked 
image and its mask image. Most of the stacked image has 
two exposures, an overlapping area has more than three ex¬ 
posures, and there are some small holes located in the CCD 
gaps that are blank. Some of the CCD gaps have only one 
exposure, so the depth is about 0.75 mag shallower in those 
locations. 

8. Photometry 

The SCUSS photometric pipeline generates comprehen¬ 
sive catalogs using multiple photometric techniques so that 
different users can choose the one that suits their needs. 
Aperture, automatic aperture, point-spread function (PSF), 
and model photometry are applied to both SCUSS stacked 
and single-epoch images. The model photometry is consis¬ 
tent with the SDSS model photometry that utilizes the SDSS 
r-band model shape parameters to measure brightnesses on 
the SCUSS u-band images. 

8.1. Source detection 

Source detections for astronomical science images are 
never complete, especially at the fainter magnitude end. 
Morphologies of extended sources in u band are more frag¬ 
mented and diffuse than those in other optical bands. There¬ 
fore, many fainter sources with low surface brightness might 
be missing. In addition, most of the science projects based 
on SCUSS data also need to use deeper and redder other 
bands from other large-scale surveys. More than 3/4 of the 
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Fig. 12.— Histograms of seeing (in arcsec), airmass, photometric zeropoint (in mag for 1 ADU s 
mag arcsec"^) for all SCUSS u-band observations. 
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Fig. 13.— Examples of a stacked image (left) and its mask image (right). The central coordinate is a = 44°.7536, i5 = 6°.1074. 
North is up and east is left. A galactic cluster of Abell 400 (marked as a plus symbol in the left plot) happens to appear in this 
area. The mask image shows the exposure numbers. There are two black boxes presenting no observations. Most of the area 
has two exposures, a small part shows more than two exposures, and some CCD-gap areas have only one exposure. The small 
horizontal lines are caused by the missing CCD rows that were not recorded by the camera controller during 2011 and 2012. 
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SCUSS area is covered by the SDSS. Our sources include 
both SDSS detected objects and detections unique to our 
SCUSS images. For the area not covered by the SDSS, we 
perform photometry for only SCUSS detections. The result¬ 
ing catalog within these areas will be useful to match with 
other wide imaging surveys (e.g. Pan-STARRS). 

The SDSS objects with any one of the PSF, Petrosian, 
model, and CModel ugriz-hand magnitudes brighter than 
23.5 mag are selected. SExtractor is used to detect objects 
in SCUSS stacked images, and these are matched with the 
SDSS objects. The mismatched objects are rematched with 
SDSS catalogs with magnitudes fainter than 23.5 mag in 
order to hnd missing SDSS objects. The rest of the objects 
are SCUSS unique detections. 

8.2. Aperture photometry by DAOPHOT 

Circular aperture photometry is a simple procedure to 
measure magnitudes of sources. The core code of the aper¬ 
ture photometry in DAOPHOT is utilized to measure aper¬ 
ture magnitudes (Stetson 1987). We use 12 apertures with 
radii ranging from 1".4 to 18" (from 3 to 40 pixels; see Table 

4). 

Table 4: Aperture radii for the aperture photometry 


No. 

Radius in Pixels 

Radius in Arcsec 

1 

3 

1.36 

2 

4 

1.82 

3 

5 

2.27 

4 

6 

2.72 

5 

8 

3.63 

6 

10 

4.54 

7 

13 

5.90 

8 

16 

7.26 

9 

20 

9.08 

10 

25 

11.35 

11 

30 

13.62 

12 

40 

18.17 


It is important to apply aperture corrections to aperture 
magnitudes due to flux loss within a finite aperture size. 
For smaller apertures and larger seeing, this is especially 
important. Figure 14 presents growth curves under differ¬ 
ent seeing conditions. The growth curve is calculated as the 
magnitude difference between one of 12 apertures and the 
7".3 radius aperture as function of aperture radius. This 
reference aperture is the one used for photometric flux cali¬ 


bration. We choose isolated and point-like objects with pho¬ 
tometric errors in all apertures less than 0.05 mag to derive 
the growth curve. Outlier objects are eliminated when the 
median growth curve is calculated. 

8.3. Automatic photometry by SExtractor 

SExtractor (Bertin & Arnouts 1996) provides precise 
magnitudes of sources using automatic aperture photome¬ 
try, generating so-called ” automatic magnitude”, which is 
motivated by the Kron algorithm (Kron 1980) An ellipti¬ 
cal aperture is automatically determined for each object to 
integrate the flux. Although most of the flux is expected 
to lie within the elliptical aperture and the flux loss should 
be almost independent of the source magnitude, we discover 
that in fact, the flux loss does change with both the source 
magnitude and seeing. The magnitudes of objects that are 
brighter or objects observed with worse seeing have more 
corrections. We consider an equivalent circle whose area is 
equal to that of the ellipse. Thus, the radius of this cir¬ 
cular aperture is equal to the root of the product of the 
semi-major and semi-minor axis lengths. We call this cir¬ 
cular radius the characteristic radius, which has the same 
meaning as the aperture size in normal aperture photome¬ 
try discussed above. Figure 15 presents the aperture cor¬ 
rections for the automatic magnitudes with different seeing. 
The black dots are point-like and isolated objects detected 
by SExtractor. They lie along with the growth curve of the 
aperture correction described in the previous section. Thus, 
an interpolation from the curve is good enough for correcting 
automatic magnitudes. 

We also perform Petrosian-like photometry on SCUSS im¬ 
ages by SExtractor. The Petrosian aperture is very similar 
to the Kron one. They share the same position angle and 
ellipticity. The Petrosian aperture radius is determined by 
the ratio of the isophotal brightness at a certain radius and 
average surface brightness within this radius. The ratio is 
set to be 0.2 and the corresponding Petrosian radius is larger 
than the Kron radius. The aperture corrections for Petrosian 
magnitudes can be estimated by the same way as mentioned 
above. 


^ https: / / WWW. astromatic. net / software / sextractor 
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Fig. 14.— Growth curves for aperture corrections at different seeing. The black pluses are growth curves for about one hundred 
isolated and point-like objects. The red lines with error bars are the medians and standard deviations of the corrections for 
different apertures. The last panel shows the median growth curves in the previous three panels in order to present the impact 
of seeing. 
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Fig. 15.— Difference between the automatic and circular aperture magnitudes of 7".i radius as a function of the characteristic 
radius. Black dots are isolated point sources. Growth curves for corresponding circular aperture corrections shown in Figure 14 
are overlaid in red. The cyan curves are the cubic spline interpolation of those growth curves. 
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8.4. PSF photometry 

The PSF is obtained using PSFEx^ (Bertin 2011). The 
form of the PSF in PSFEx is expressed as a linear combina¬ 
tion of basis vectors. The pixel basis is selected in our PSF 
modeling. The spatial PSF variation on the focal plane usu¬ 
ally shows a smooth profile, which can be modeled by a low- 
order polynomial. For SCUSS single exposures, a second- 
order polynomial is good enough to describe the PSF vari¬ 
ation over the CCD plane. For SCUSS stacked images, a 
seventh-order polynomial is used to describe the complexity 
of the image quality. 

A code is specially designed to measure the PSF mag¬ 
nitudes at known object positions. The code takes the 
SCUSS image, SDSS and SCUSS-only object positions, the 
bad-pixel list, and the PSF model derived by PSFEx as in¬ 
puts and then outputs the Gaussian-fitted position, local 
sky background and its error, local PSF FWHM, fitted PSF 
integrated flux and its estimated error, and a flag tagging 
the status of each object. Figure 16 shows the flowchart 
of the PSF photometry code. The local sky background for 
each object is measured using the pixels at r > 7.5 xFWHM. 
Here, the FWHM is the full width at half maximum of the 
local PSF model. Outliers, such as cosmic rays and sig¬ 
nals from real objects, are iteratively rejected by a sigma- 
clipped algorithm. After the sky background is subtracted, 
the object position on the CCD is fitted by a two-dimensional 
Gaussian function. The position is allowed to shift because 
the coordinates in SDSS redder bands are slightly different 
from those in the SCUSS u band due to atmospheric refrac¬ 
tion and star proper motions. The pixels centered at the 
initial position with r < 2.5xFWHM are considered in the 
calculation. If the new fitted position is more than 2.0 pixels 
away from the old one, the following PSF photometry will be 
performed at the original position. The PSF model is inter¬ 
polated to the same pixel scale of the CCD image and fitted 
to the fluxes of the object pixels within r = 2.Ox FWHM. 
Flags are also provided by our PSF photometry to show the 
reliability. They are coded in decimal and expressed as a 
sum of powers of 2: (1) CCD artifacts; (2) bad pixels; (4) 
including saturated pixels; (8) contaminated by neighbors; 
(16) near image edges. 

We compare our PSF measurements with one of the popu¬ 
lar PSF photometry software packages (DAOPHOT), which 


® http: / / WWW. astromatic. net / software / psfex 


is a widely used package for accurate stellar photometry de¬ 
signed to deal with crowded fields. We use DAOPHOT to 
find objects in one single-epoch image and at the same time 
give the PSF magnitude measurements. Then, our code per¬ 
forms PSF photometry for those objects with the PSF model 
derived by PSFEx. The magnitude comparison is shown in 
Figure 17. Two groups of objects are chosen: a brighter 
one with 16 < u < 18 mag (left in Figure 17) and a fainter 
one with 21 < u < 22 mag (right in Figure 17). The scat¬ 
ter of the PSF magnitude differences between our code and 
DAOPHOT are 0.009 mag for the bright group and 0.021 for 
the faint group, respectively. Our PSF photometry is con¬ 
sistent with DAOPHOT, although the PSFs used by these 
two techniques are derived in different ways. 

8.5. Model photometry 

The SDSS also provides a type of model magnitude, 
or ’’modelMag,” which can measure the unbiased colors of 
galaxies in the absence of color gradients through equivalent 
apertures in different bands. The ’’modelMag” is generated 
by choosing a shape model, either a deVaucouleurs or expo¬ 
nential profile, based on its best-fit likelihood in the SDSS 
r band, then convolving them with seeing in other bands, 
and finally forcing magnitude measurements with the same 
aperture shape. 

We divided SDSS objects into two groups, point sources 
and extended sources, based on the SDSS classification. For 
extended sources, we construct their theoretical 2D models 
with the effective radii, axis ratios, and position angles from 
the SDSS measurements. These models are convolved with 
local PSF profile derived by PSFEx. For point sources, we 
use local PSF from PSFEx directly as their models. In ad¬ 
dition, the extended sources with small sizes or low bright¬ 
nesses (effective radii less than 0.5 arcsec or SDSS r-band 
magnitude fainter than 23.5 mag) are also treated as point 
sources. 

After the models are constructed, model amplitudes are 
calculated as the ratios of the models and raw SCUSS fluxes. 
The pixels within r = 1.0 x FWHM are considered in the 
calculation, which is optimized for faint SDSS extended 
sources. According to the amplitudes, we then compute 
the corresponding deVaucouleurs and exponential magni¬ 
tudes. Therefore, the ’’modelMag” in SCUSS u band can 
be estimated by the SDSS r-band deVaucouleurs and expo¬ 
nential profiles. The SCUSS model magnitude is aperture- 
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Fig. 16.— Flow chart of our PSF photometry pipeline. Firstly, the PSF model is extracted by combining SExtractor and 
PSFEx. Secondly, the coordinate for each object is fitted by a Gaussian function. Thirdly, the image of each object with the 
sky background subtracted is fitted by the local PSE profile. Finally, the photometric results together with the fitted positions 
and flags are given as outputs. 
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Fig. 17.— Comparisons of our PSF photometry with DAOPHOT. Left is the comparison for common bright objects with 
16 < M < 18 and right is the one for common faint objects with 21 < u < 22. 
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corrected to make the magnitudes of bright point sources 
(16 < u < 20) equal to the SDSS w-band model magnitudes. 

8.6. Co-added photometry 

The stacked images are composed of single-epoch images 
taken with different observational conditions. The seeing 
varies between each single-epoch image, which makes the 
PSF profile of the stacked image quite complicated. Sub¬ 
sequently, the fraction of flux loss due to finite apertures 
for objects in different parts of the stack images might be 
different. In addition, the PSF profile cannot be perfectly 
determined unless the seeing of related single-epoch images 
is similar enough. On the other hand, the PSF profile of a 
single image varies smoothly and it is much easier to model. 
We can first perform photometry for an object observed with 
different exposures and then co-add its fluxes to generate dif¬ 
ferent co-added magnitudes at the catalog level. 

The aperture photometry for single-epoch images is ob¬ 
tained by DAOPHOT with 12 apertures as defined before. 
Aperture corrections are applied to the resulting aperture 
magnitudes. Corresponding corrected fluxes are weighted by 
the errors and averaged to generate the co-added aperture 
magnitudes. 

The PSF photometry for single-epoch images is performed 
using our PSF fitting code. The PSF magnitudes are also 
aperture-corrected by comparing with the aperture magni¬ 
tudes in a 7".3 radius. The co-added PSF magnitudes are 
calculated by averaging the fluxes weighted by their errors. 

We perform the model photometry for single-epoch im¬ 
ages and measure the exponential and deVaucouleurs mag¬ 
nitudes. We adopted model magnitude with a higher SDSS 
r-band likelihood. Aperture corrections are applied to make 
model magnitudes equal to the SDSS model magnitudes in 
the case of unresolved objects. The co-added model magni¬ 
tude is calculated by averaging the model fluxes weighted by 
their errors. 

9. Some photometric comparisons 

The photometric methods described above are applied to 
stacked and single-epoch images. Figure 18 illustrates the 
general photometric scheme. The detections are based on 
both SDSS and SCUSS images. We perform aperture, PSF, 
model, and automatic photometry for the stacked images. 
Note that the automatic magnitudes are only based on ob¬ 


jects solely detected by SExtractor. We also obtain co-added 
aperture, PSF, and model magnitudes from photometry of 
single-epoch images. Since point sources can be best mod¬ 
eled while extended sources are much more complicated, the 
photometry of point-like objects is better suited for compar¬ 
isons of our different photometric methods. Thus, the com¬ 
parisons of point sources are presented for most cases below 
and, for comparison with SDSS, we correct the magnitude 
with the SCUSS/SDSS color term. 

9.1. Photometry of stacked images 

We compare the different SCUSS photometric measure¬ 
ments for stacked images with the SDSS PSF or model mag¬ 
nitudes of point sources in Figure 19. The SCUSS model 
magnitude is compared with the SDSS u-band model mag¬ 
nitude, while other magnitudes are compared with the SDSS 
PSF magnitude. We choose the 5 pixel (2".27) to present 
the aperture magnitude in this figure. For the smallest aper¬ 
tures, the aperture magnitudes might be problematic due to 
inhomogeneous image quality in some stacked images. Table 
5 gives the magnitude offset and scatter in two magnitude 
intervals: 16 < it < 19 mag and 20.5 < m < 21.5 mag (i.e. 
^21 mag). The automatic magnitude is best among all pho¬ 
tometric magnitudes for stacked images, with a scatter of 
0.033 for bright stars and 0.178 for faint stars at it ~ 21.0 
mag. 

9.2. Co-added photometry for single images 

The comparisons of the co-added PSF, aperture, and 
model magnitudes with the SDSS PSF or model magnitudes 
are shown in the right panels of Figure 19 and the last three 
rows of Table 5. The co-added PSF magnitude performs 
the best for point sources among all the magnitude types. 
The scatter is 0.033 for brighter sources and 0.174 for fainter 
sources at it ~ 21.0 mag. The co-added aperture magnitude 
is also adequate if a proper aperture radius is considered. 
The aperture size should be chosen according to the object 
type, object brightness, signal-to-noise requirement, etc. 

9.3. Comparisons with the CFHTLS deep it band 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS; Astier et al. 2006) used the wide-held optical 
imaging camera on CFHT to obtain deep multicolor pho¬ 
tometry over wide areas over 5 years.The photometric sys¬ 
tem is similar to that of the SDSS, except the it Hlter. As 
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Fig. 18.— SCUSS photometric scheme. The output magnitudes include the automatic, Petrosian, aperture, and PSF magnitudes 
for the stacked images and co-added aperture and PSF magnitudes for the single-epoch images. 
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Fig. 19.— Magnitude differences between different SCUSS magnitudes and the SDSS PSF or model magnitudes of points sources 
as functions of SCUSS magnitudes. The first-row panels and the rest of the left panels are comparisons of automatic, Petrosian, 
PSF, and aperture magnitudes for the stacked images. The right panels from the second row to bottom are comparisons of 
co-added PSF, aperture, and model magnitudes for the single-epoch images. The aperture magnitude shown here is the one with 
an aperture radius of 4".5. The two lines above and below zero show ±1 ct scatter around the average. 
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Table 5: Magnitude offset and scatter of point sources for the 
SCUSS magnitudes compared with the SDSS PSF or model 
magnitudes 


Magnitude 

16 < 

M < 18 

20.5 < 

u < 21.5 

Offset 

Scatter 

Offset 

Scatter 

Automatic (stacked) 

0.001 

0.033 

0.008 

0.178 

Petrosian (stacked) 

0.002 

0.033 

0.065 

0.200 

PSF (stacked) 

-0.003 

0.047 

0.025 

0.175 

Aperture (stacked) 

-0.002 

0.044 

0.017 

0.176 

Model (stacked) 

-0.006 

0.059 

0.004 

0.177 

Co-added PSF 

-0.004 

0.033 

0.021 

0.174 

Co-added aperture 

0.003 

0.034 

0.010 

0.177 

Co-added model 

-0.003 

0.049 

0.003 

0.175 


SDSS scatter for brighter sources is similar to the SCUSS 
scatter, while it is much larger than the SCUSS scatter at 
the faint magnitude ends. The histograms in the right pan¬ 
els of Figure 20 show the magnitude difference distributions 
with 22 < ucFHT < 23. The scatter for the SCUSS and 
SDSS are also presented in these panels. The photometric 
accuracy of the SCUSS is much better than that of the SDSS 
for fainter sources due to deeper imaging. 

9.4. Photometry for extended sources 

Extended sources are much more complicated, especially 
in ultraviolet bands, since their morphologies look much 
more fragmented than in redder bands. By comparing dif¬ 
ferent magnitude measurements on stacked images with the 
SDSS model magnitude, we calculate a scatter at u ~ 21 
of about 0.2, 0.31, 0.4, 0.23 for automatic, aperture (2".27), 
PSF, and model magnitudes, respectively. The co-added 
photometric magnitudes give similar results. The automatic 
magnitudes seems to be the best for extended sources, but 
they only describe the brightnesses of SCUSS-detected ob¬ 
jects. The model and aperture magnitudes are also ade¬ 
quate. The PSF photometry fails to measure the magnitudes 
of galaxies. 

9.5. Guidelines to use magnitudes 

The choice of photometric technique is dependent on the 
science one wants to do. Here we present some general guide¬ 
lines. More than 90% of the stacked images are assembled 
from single-epoch images with consistent image quality. The 
photometry on these images is as good as the corresponding 
co-added photometry. However, as the background noise in 
single-epoch images is larger than that in stacked images, 
there are about 23% more objects with available magni¬ 
tude measurements in stacked images than with available 
co-added magnitudes. 

The automatic magnitude performs as well as the PSF 
and model magnitudes because it can adaptively fit elliptical 
apertures to both point and extended sources with similar 
flux loss. It can be regarded as a universal magnitude for 
both point-like and extended sources. But unlike other pho¬ 
tometric methods, the automatic photometry is based on the 
objects detected on SCUSS images, which are about 35% of 
SDSS objects with available SCUSS w-band fixed-parameter 
measurements. For point sources like quasars and stars, 
. ~ i , , , the PSF and aperture magnitudes with appropriate aper- 

^http: //cmt. hawaii.edu/Instruments/Imaging/MegaPrime/generalinformation.html 


one component of the CFHTLS, the wide synoptic survey 
covers 155 square degrees in four patches down to i =24.5. 
The rt-band magnitude limit in the wide survey reaches 25.2 
mag, which is much deeper than the SCUSS u band. Thus, 
the catalogs of the wide survey can be regarded as a refer¬ 
ence when comparing the SCUSS photometry with that of 
the SDSS. There are two CFHTLS wide regions, W1 and 
W4, overlapping with the SCUSS footprints. The W4 field 
(a = 22‘'13“18%5 = -t01‘^19“00") has an area of 25 deg^ 
fully covered by the SCUSS and it is located at a higher 
declination, where the SCUSS data quality is more typi¬ 
cal than it is in the W1 field. We make some photometric 
comparisons by using the W4 catalogs. 

Figure 20 shows the photometric comparisons of point 
sources between the SCUSS and SDSS with the CFHTLS 
wide data as a reference. The left panels in this figure 
present the magnitude difference between the SDSS (in blue 
points) or SCUSS u-band magnitude (in red points) and the 
CFHTLS automatic magnitude as a function of the CFHTLS 
u magnitude. The SDSS magnitude is the PSF magnitude 
and the SCUSS magnitudes are automatic, co-added PSF, 
and co-added aperture (2".27) magnitudes from top to bot¬ 
tom, respectively. These three types of SCUSS maguitudes 
are considered to be the best flux measurements for point 
sources. The number of objects for the automatic magni¬ 
tude is less than those of the other two magnitude types 
(mainly at the faint end) due to different source detection by 
SExtractor itself. Both the SCUSS and SDSS it-band mag¬ 
nitudes are converted to the CFHTLS photometric system 
by the color term as indicated in the CFHTLS webpage 
MCFHT = MsdSS/SCUSS “ 0-241('«SDSS/SCUSS “ SSDSs)- The 
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Fig. 20.— Photometric comparisons of point sources between the SCUSS (in red) and SDSS (in blue) with the CFHTLS W4 
data as a reference. The SDSS magnitude is the PSF magnitude and the CFHTLS magnitude is the automatic magnitude. The 
left panels are the magnitude differences between the SCUSS or SDSS and the CFHTLS as a function of the CFHTLS u-band 
magnitude, and the right ones are the normalized distributions of the differences at the faint magnitude end (the hatch area: 
22 < u < 23 mag). From top to bottom, they are comparisons of three types of SCUSS magnitudes: the automatic, co-added 
PSF, and co-added aperture magnitudes, respectively. Only about 20,000 objects are randomly selected in this figure in order 
to avoid crowding. The color texts in the right panels gives the magnitude scatter of the SCUSS and SDSS, which corresponds 
to the 68.3% confidence level. 
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ture sizes are recommended. The aperture size needs to be 
determined based on scientific objectives. For nearby galax¬ 
ies with extended morphological structures, the automatic 
magnitude and aperture magnitude are good choices. The 
SCUSS model magnitude is defined the same as the SDSS 
” modelMag”. When combining magnitudes of other SDSS 
bands to measure the colors of extended sources, it is better 
to use the model magnitude. 

10. Summary 

The SCUSS survey is a wide-field u-band sky survey in 
the Southern Galactic cap. The survey used the Bok tele¬ 
scope on Kitt Peak and the filter is close to the SDSS u band. 
The survey observations were completed by the end of 2013 
and the total area is about 5000 deg^. This paper describes 
the detailed data reduction dedicated to the survey, includ¬ 
ing basic image processing which has some special features 
related to the detectors, astrometric and photometric cali¬ 
brations, and photometry. The general astrometric error is 
about 0".13. The SCUSS photometric calibration is tied to 
the SDSS catalogs and is performed for each amplifier of a 
CCD due to gain and weather variations. 

We apply different photometric techniques to the stacked 
images, including automatic photometry by SExtractor, 
aperture photometry by DAOPHOT, PSF photometry, and 
model photometry. Our PSF photometry with more con¬ 
trollable parameters is consistent with DAOPHOT. The 
model photometry is similar to the SDSS ’’ModelMag”, 
which uses the SDSS r-band model-derived shape param¬ 
eters and SCUSS PSF profiles to make consistent and unbi¬ 
ased model magnitude measurements. We perform photom¬ 
etry on stacked images and also on single-epoch images, from 
which co-added photometry is derived. More than 90% of 
stacked images are assembled from single-epoch images with 
consistent quality. There are about 23% more objects with 
available magnitudes on stacked images than with available 
co-added magnitudes. The photometry on these stacked im¬ 
ages is as good as the co-added photometry. However, for 
the rest of the stacked images, their photometry is worse 
than the co-added photometry due to uneven image quality. 
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